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ABSTRACT 

In this paper we present results of an analysis of radio-quiet QSO environments. The 
aim is to determine the relation between QSOs, galaxies and the mass distribution as 
a function of redshift. We cross-correlate a sample of ~ 150 QSOs from optically and 
X-ray selected catalogues with faint, bj < 23, galaxies. These data allow us to probe 
the galaxy clustering environment of QSOs out to z ~ 1.0 — 1.5. Far from giving a 
positive correlation, at z < 1.5 the QSO-galaxy cross-correlation function is marginally 
negative, with oj(6 < 120") = —0.027 ± 0.020. Colour information suggests that the 
anti-correlation is most significant between the QSOs and the red galaxy population. 

We have constructed models to predict the QSO-galaxy cross-correlation, using 
the known form of the galaxy N(z) at bj < 23, and assuming a variety of clustering 
evolution rates. Cases in which QSOs exist in rich cluster environments are comfortably 
ruled out at more than 5<r and the results are more consistent with a 'normal' galaxy 
environment for radio-quiet QSOs. If the small anti-correlation is interpreted as an 
effect of gravitational lensing, this conclusion is not altered. In this case, the data 
are only ~ la below the low clustering amplitude models, while the high amplitude 
models are still comfortably rejected. We therefore conclude that these QSOs may not 
be much more highly biased than optically selected galaxies. 

Key words: galaxies: clustering - galaxies: evolution - quasars: general 



1 INTRODUCTION 

Relating QSOs to the underlying mass distribution and to 
other mass tracers is a vital step in the process of using QSOs 
to measure large-scale structure in the Universe. In order 
to determine cosmological quantities, such as the density 
parameter flo, from the evolution of large-scale structure we 
are required to know the evolution of the mass distribution. 
Therefore, knowledge of the biasing of QSOs with respect to 
this mass distribution is of the utmost importance, and the 
relationship between galaxies and QSOs is an essential step 
in understanding this bias. The local environments of QSOs 
also give important clues to the physical processes which 
form these objects and control their evolution. 

Direct imaging of the regions around QSOs is the most 
obvious way to investigate their environments. Yee & Green 
(1984) imaged 3' x 3' fields centred on individual QSOs over a 
wide range in redshift, (0.05 < z < 2.05). They found signif- 
icant numbers of excess galaxies (to a limiting magnitude of 
~ 21 in the Gunn r band) associated with QSOs at z < 0.5. 
The luminosity distribution of the excess galaxies strongly 
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suggested that these QSOs are associated with galaxies at 
the distance i mplied by their cosm ological redshift. Further 
investigations ( Yee fc Green 1987 ) showed that at z ~ 0.6 
some radio-loud QSOs were found in environments as rich 
as those of Abell class 1 clusters. Significant evolution was 
also detected from z = 0.4 to z — 0.6 in this radio-loud 
sample, with the QSO-galaxy covariance amplitude increas- 
ing by a fact or of ~ 3 as redshi ft increases. An extension 



of this work (Ellingson et al. 1991) found that radio-quiet 
QSOs by contrast exist in significantly poorer environments 
than their radio-loud counter-parts. The results concerning 
radio-loud QSOs have been confirmed by several other in- 
dependent studies (e.g. Hintzen et al. 1991) which found 
significant numbers of excess galaxies associated with QSOs 
out to z ~ 1.5. However, observations of radio-quiet QSOs 
have found a range of results generally consistent with these 
objects existing in poorer environments, similar to that of 
an average galaxy. Smith, Boyle & Maddox (1995) find that 
the cross-correlation function between low redshift (z < 0.3) 
X-ray selected QSOs and 6,j < 20.5 galaxies from the APM 
Galaxy Survey (Maddox et al. 1990) is consistent with the 
auto-correlation function of the APM Galaxy Survey, sug- 
gesting that this population of QSOs is unbiased with re- 
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Table 1. Details of the AAT photographic plates used in our 
QSO-galaxy cross-correlation 



spect to galaxies. Extending this work, Smith (1998) finds 
that X-ray selected QSOs to z = 0.7 exist in poor environ- 
ments. Boyle & Couch (1993) have found no excess galaxy 
population associated with radio-quiet QSOs at z ~ 1, while 
Hutchings, Crampton & Johnson (1995) find that a number 
(9 out of 14) of radio-quiet QSOs exist in compact groups of 
star forming galaxies at z = 1.1. Clearly there is a difference 
between radio-loud and radio-quiet QSOs, with radio-loud 
QSOs inhabiting richer environments. Radio-quiet QSOs ap- 
pear to inhabit environments similar to those of normal 
galaxies, although the Hutchings et al. result is somewhat 
discrepant from this hypothesis. In this paper we will fur- 
ther investigate the problem of QSO environments with a 
large (~ 150) sample of optically and X-ray selected QSOs 
covering a wide range in redshift. The second paper in the se- 
ries (Croom & Shanks 1998, in preparation) will use deeper 
imaging to constrain QSO environments at higher redshift. 

In Section 2 we describe our data and the methods em- 
ployed. We present our cross-correlation results in Section 
3, and compare them to model predictions in Section 4. A 
discussion of our results and our conclusions are contained 
in Sections 5 and 6. 



2 D ATA AND METHODS 



The faint galaxy samples were taken from deep AAT plates 
in five 40' x 40' fields. Details of the plates are given in Ta- 
ble |l|. Each field has a fcj plate and two fields also have tf 
plates. The plates were scanned by the COSMOS measuring 
machine. Details of the analysis of the first 3 fields (SGP2, 
SGP4, QSF3), including star-galaxy separation and photo- 
metric calibration are given in Jones et al. (1991). Fields 



F855 and F864 (Roche et al. 1995; Roche 1994) were simi- 



larly analysed. The plate scale for all five fields is 15.2 arc- 
sec/mm. The magnitude limit for the bj plates is ~ 24 mag 
and the completeness limit is bj — 23.5. We use the galaxies 
to a limiting magnitude of bj = 23.0, which is comfortably 
brighter than the completeness limit. At this magnitude the 
errors in the star-galaxy separation are ~ 5%. 

The QSOs in this analysis are from both optically and 
X-ray selected samples. Fields SGP2, SPG4 and QSF3 con- 
tain QSOs from the Durham/ A AT UVX selected sample 
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Figure 1. The redshift distribution of all the QSOs, optically and 
X-ray selected, used in this analysis. The total number of QSOs 
is 152. 



QBoyle et al. 1990Q , limited to B = 21 mag. The F855 and 
F864 fields conta in QSOs from the deeper ESQ /AAT colour 
selected sample (Boyle et al. 1991; Zitelli et al. 199S ) with a 



limiting magnitude of B — 22.5. The X-ray selected QSOs 
were taken from optical follow-up o f a deep ROSAT survey 
dAlmaini 199fj ; |Shanks et al. 1998] ). The QSOs are > 4<r 
detections with the R OSAT Position Sensitive Proportional 
Counter (PSPC). Details of the QSO samples used are given 
in Table [| and the redshift distribution of the combined 
sample of 152 QSOs is shown in Fig. [j]. It should be noted 
that a large number of QSOs were selected independently 
using both optical and X-ray techniques. As pointed out in 
Section 1, radio-loud and radio-quiet QSOs appear to have 
different environments. We have therefore attempted to de- 
termine whether any of our optically and X-ray selected 
QSOs are, in fac t, radio-loud. Using the NRAO VLA Sky 
Survey (NVSS) ( |Condon et al. 1998| ) at 8 > -40° and the 
Parkes-MI T-NRAO (PMN) Survey at 5 < -40° (|Gregory 



et al. 1994) we find that none of our QSOs are within 30 of 
a radio source at the ~ 2.5mJy level (NVSS) and ~ 25mJy 
(PMN). This is somewhat surprising, as we expect ~ 5% of 
optically selected QSOs to be radio-loud. 

The QSO coordinates (optical counterpart positions in 
the case of the X-ray selected objects) were transformed 
to the x,y system of the COSMOS scanned fej plates using 
the Starlink ASTROM software; ~ 40 stars were used to 
produce the 6 coefficient transform. The RMS residuals for 
these transforms are shown in Table ^, all are less than 1". 
The five fields are plotted in Fig ^, showing all galaxies to 
fej = 23 with the positions of the QSOs overlaid on top of 
this. 

The cross-correlation was carried out in the standard 
manner, with the galaxy distribution being compared to a 
random distribution of points with the same angular selec- 
tion function (that is avoiding holes) as the galaxies. The 
number density of randoms was 100 times the density of 
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Figure 2. The x,y positions of galaxies to ft, j = 23 from COSMOS scans of AAT photographic plates, together with the positions of 
the QSOs. a) Field SGP2, b) field SGP4, c) field QSF3, d) field F855, e) field F864. Open squares arc the optically selected QSOs while 
crosses are the X-ray selected QSOs, The regions removed due to plate defects and bright stars can be seen. It should be noted that a 
number of QSOs are both optically and x-ray selected. 



galaxies. The angular cross-correlation function, w qg (0), is 
defined as 



_ N qe (9)N r 



(1) 



where N T is the total number of random points, N g is the 
total number of galaxies, and iV q g and iV qr are the num- 
ber of QSO-galaxy and QSO-random pairs with separation 
9 ±59/2 respectively. This was calculated in concentric an- 



nuli of width 30". The inner 5" is removed to avoid any con- 
fusion between very near companions and the QSOs them- 
selves, and any QSOs which might have been classified as 
galaxies on the COSMOS scans. The errors shown in the 
figures below are Poisson errors: 



AW : 



w(0) + 1 



(2) 



To test the Poisson error estimates we carry out two tests; 
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Figure 2. Cont. 



Table 2. Details of the QSO samples used in this analysis. It 
should be noted that many QSOs were selected by both optical 
and X-ray methods. 



Field 
Name 


No. of 
Opt. QSOs 


No. of 
X-ray QSOs 


Total No. 
of QSOs 


r.m.s. 
Ast. (") 


SGP2 


20 


20 


30 


0.86 


SGP4 


10 


27 


33 


0.43 


QSF3 


17 


21 


26 


0.39 


F855 


21 


19 


31 


0.26 


F864 


22 


20 


32 


0.69 



the first is to replace the QSOs with a random selection of 
galaxies and carry out a cross-correlation between these ran- 
dom galaxies and the rest of the galaxy sample. A variance is 
calculated from 100 realizations of this process. The second 
test calculates the cross-correlation between the QSOs and 
a bootstrap re-sampled galaxy distribution, again using 100 
realizations to find a variance. For all QSOs with z < 1.5 the 
Poisson error on the amplitude within 120" is ±0.020 (see 
Section |^ below) . The error estimated from replacing QSOs 
with random galaxies is ±0.026 while the bootstrap estima- 
tor gives ±0.020. The excess variance from the first test is 
due to the positive galaxy auto-correlation. The bootstrap 
method uses the QSO-galaxy cross-correlation which, as we 
show below, has no significant clustering signal. 

The background density of galaxies is determined from 
each field separately which forces the integral of uj(9) to be 
zero at the largest scales studied. This is the well known 
'integral constraint', which is simply due to the fact that 
there is still appreciable clustering of galaxies on the scale 
of our fields. In our case we do not know a priori that there 



will be a significant signal on scales of ~ 20', therefore an 
accurate estimate of the integral constraint is not possible. 
Given this, we can still determine a possible integral con- 
straint correction in order to investigate its possible effect 
on our results. The angular correlation function determined 
in an area f2 is biased low by the amount 



1 



w(6>i 2 )dftidfi 2 



(3) 



We then assume the standard power-law for the cross- 
correlation function, u>(9) = A9~ 08 . We find the following 
integral constraint for the five fields: 3.29 A (SGP2), 3.42,4 
(SGP4), 3.364 (QSF3), 3.044 (F855) and 2.894 (F864). 
The auto-correlation function amplitude for 6j ~ 23 galax- 
ies is ~ 0.002. If we were to assume that this value was 
a reasonable estimate of the QSO-galaxy cross-correlation 
then the integral constraint in our fields would be ~ 0.006. 
We will demonstrate below that the addition of an integral 
constraint of this amplitude has no significant effect on our 
conclusions. 



3 CROSS-CORRELATION RESULTS 

We now present the results of the cross-correlation anal- 
ysis between both the optically and X-ray selected QSOs 
and the galaxy catalogues. First, we look at the angular 
cross-correlation function for each field individually; these 
are shown in Fig. ^a-e while Fig. ||f shows a combination of 
all five fields. We show the results for all the known QSOs 
within each field (X-ray and optically selected QSOs com- 
bined) . From these plots it is clear that there is a significant 
anti-correlation in the F855 and F864 fields while the SGP4 
field is the only one which appears to show a positive cor- 
relation. The combined cross-correlation function therefore 
shows a small negative signal on small scales. Fig. ^| was de- 
rived using all the QSOs in each field (including the whole 
range of redshifts and both X-ray and optically selected ob- 
jects). In order to look in detail at the source of the signal 
we calculate below a clustering amplitude for each individual 
QSO. 

We calculate the individual cross-clustering amplitude 
for each QSO in a non-parametric fashion, simply using the 
integrated excess of galaxies out to a defined radius which we 
choose to be 120", giving u>(9 < 120") (still using an inner 
limit of 5"). A radius of 120" is equivalent to 0.5 /i -1 Mpc 
in proper co-ordinates at a redshift of z ~ 1 for an Qo — 1 
flat Universe. We can relate this measure of clustering to a 
parametric clustering amplitude assuming the simple power- 
law form for the angular cross-correlation function, 



The number of QSO-galaxy pairs between 9\ and #2 is 



< e < . 



2N S 



02 



(4) 



(•>) 



where N E is the number of galaxies. This can then be used 
to relate the non-parametric ui{9 < 120") to the clustering 
amplitude A qs via Eq. [I] for 5" < 9 < 120". We assume a 
power- law slope with 7 = 1.8. 

In Fig. [| we show the clustering amplitudes for individ- 
ual QSOs within 120" as a function of redshift for each field. 
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Figure 3. The QSO-galaxy angular cross-correlation function for each of the five AAT fields individually, including all QSOs. f) shows 
a pair weighted combination of all five fields. 
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Figure 4. The QSO-galaxy clustering amplitude for individual QSOs in each AAT field (both optically and X-ray selected QSOs) as a 
function of rcdshift. f) shows all fields combined. Individual clustering amplitudes are measured within a radius of 120". 
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Figure 5. QSO-galaxy clustering amplitudes within 120" for QSOs as a function of redshift. a) individual X-ray selected QSOs and b) 
individual optically selected QSOs. c) and d) show the same data binned as a function of redshift with Az = 0.25. 



These show that the positive correlation in the SGP4 field is 
mainly due to three QSOs with strong positive signals. One 
of these is at z = 2.861 and has w(0 < 120") = 0.62 ± 0.20 
(equivalent to A qg = 0.025 ±0.008). We do not expect to see 
any galaxies brighter than fej = 23 at this redshift and hence 
this is either a chance alignment or the effect of gravitational 
lensing on the QSO. Examination of this object by eye con- 
firms that there is a concentration of faint objects close to 
the QSO which could be a cluster of galaxies. The second of 
the three objects is only separated from the z = 2.861 QSO 
by ~ 1'. The third also shows a number of faint objects very 



close to the QSO. By contrast, the negative signal in fields 
F855 and F864 is seen to be due to a large number of objects 
with lo(6 < 120") ~ -0.2 to -0.5. There is no measurable 
correlation between redshift and clustering amplitude (de- 
termined from a non-parametric rank- correlation test), with 
objects over the entire range in redshift contributing to the 
negative signal. 

We should note that using galaxies from the AAT plates 
to the plate limit of fej ~ 24 does not have any effect on the 
results presented here. The clustering amplitudes are not 
significantly altered, although obviously the Poisson errors 
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are somewhat reduced by the increased number of galax- 
ies. Considering that at this magnitude limit the sample is 
incomplete and the star-galaxy separation is ill-defined, a 
limit of bj = 23 affords a more robust and reliable result. 

The clustering amplitudes for individual QSOs sepa- 
rated on the basis of their selection (X-ray or optical) are 
shown in Fig. || ft should be noted that a large fraction 
of the QSOs were selected independently by both X-ray and 
optical techniques. Of course, the X-ray selected sample also 
has an implicit optical selection imposed on it due to the 
need to identify an optical counterpart to the X-ray source 
on which to carry out follow-up spectroscopy. Fig. H shows 
that almost all of the QSOs show zero or negative correla- 
tion amplitudes, while there are a small number of outlying 
objects with significant positive signal. The only difference 
between the X-ray and optically selected QSOs is a possi- 
ble anti-correlation in the optical sample at high, z > 1.8, 
redshift. This can be seen more clearly if the results are 
combined in redshift intervals, as shown in Figs, and d. 
Here we combine the galaxy counts around each QSO in 
Az = 0.25 intervals. The X-ray selected sample (Fig. ||c) 
shows no significant signal or trend as a function of redshift. 
The optically selected sample (Fig. ||d) shows a significant 
anti-correlation at redshifts greater than z ~ f .8, with the 
other significant points (in the first and last redshift bins) 
both due to individual QSOs. At low redshift, where any 
signal due to real associations would be expected, none is 
detected; at z < 1 .5, lu{9 < 120") = -0.018 ± 0.022 for the 
X-ray selected sample and —0.051 ± 0.027 for the optical 
sample. For the combined sample of all QSOs with z < 1.5 
we obtain uj{6 < 120") = -0.027 ± 0.020 (82 QSOs). In the 
next section these results are compared to the expected val- 
ues for given clustering amplitudes and clustering evolution. 

We have colour information in two fields, with F855 
and F864 containing r-p plates. The colour distribution of 
the galaxies peaks at fej — tf ~ 1.25 so we use this limit 
to define red and blue galaxy populations. Fig. ^| shows the 
results of this analysis, with Figs. |^a and ^b showing the 
angular correlation functions for all QSOs. The red galaxies 
appear to contribute more to the anti-correlation seen in the 
F855 and F864 fields. When we look at the clustering am- 
plitudes of individual QSOs at 9 < 120" we see no obvious 
trend (Figs. []c, d), but binning these results as a function of 
redshift we see that the blue galaxies show no strong correla- 
tion or anti-correlation at any redshift, while the red galax- 
ies appear to show generally negative correlations (similar to 
that seen in the correlation of all galaxies with the optically 
selected QSOs, Fig. ^i). At z < 1.5 the correlations for the 
red and blue populations are u(6 < 120") = -0.093 ± 0.043 
and — 0.057±0.048 respectively. Further dividing the sample 
into X-ray and optically selected QSOs we see no significant 
differences, this is mostly due to the increased errors from a 
smaller number of QSOs. 

We therefore see tentative evidence for a difference in 
the cross-clustering properties of red and blue galaxies with 
QSOs. Any physical process which would cause this effect 
will be largely due to the fact that red and blue galaxies tend 
to inhabit different environments. Red galaxies are found in 
richer environments and are more strongly clustered than 
blue galaxies, which tend to inhabit the field rather than 
clusters. This hints at gravitational lensing as a mechanism 
for the anti-correlation. Below we make some estimates of 



the expected QSO-galaxy cross-correlation and relate them 
to the results found here. 



4 MODELLING QSO ENVIRONMENTS 
4.1 Limber's equation 

The angular cross-correlation function, u} a b(9), between two 
populations, a and 6, can be related to the spatia l cross- 
corre l ation function £ a b(r, z) by Limber's equation (Limber 
1953; ; [Peebles 198o| ) which can be written as 



^ab(#) 



2 J °° J °° x 4 F~ 2 Mx)M^Mr, z)dxdu 
x 2 F~ 1 ^(x)dx x 2 F- 1 ^> b (x)dx ' 



(6) 



under the assumption that the integral is dominated by ob- 
jects at almost equal cosmic distance. The F term accounts 
for varying cosmological geometries and is given by 

F=[l-(^) 2 0W)] 1/2 . (7) 

We will assume a flat cosmology with Slo = 1, that is F = 1. 
The comoving distance to a point midway between a and 
b is given by x (assuming that the separation of a and b is 
small compared to x). The small angle approximation is also 
assumed so that the proper separation between a and b is 



U i 2 2 

— rr + x a 

F 2 



1/2 



u — Xt, 



(8) 



4>n{x) and (f>b(x) are the selection functions for the two pop- 
ulations and satisfy 



(x)dx = 



N(z)dz, 



(9) 



where N(z)dz is the number of objects per unit surface area 
in the redshift interval [z,z + dz\. 

We parameterize the spatial cross-correlation function 
in the conventional power-law form: 



(10) 



U{r,z)=(j^) (l + z)-^ 3+c \ 

where r is the proper separation of two objects and ro is the 
correlation scale length at z = 0. Evolution is parameter- 
ized by e, with various values of this parameter correspond- 
ing to the following cases: e = is equivalent to constant 
clustering in proper coordinates, so called stable clustering; 
e = 7 — 3 implies clustering which is constant in comoving 
coordinates; e = 7 — 1 implies clustering which grows accord- 
ing to linear theory (for fio = 1)- We adopt the standard 
locally measured galaxy auto-correlation power-law slope of 
7 = 1.8. Although this is not necessarily correct at high 
redshift or in the case of QSO-galaxy cross-correlations, the 
measured QSO auto-cor relation is consistent with this slope 
( |Croom fc Shanks 1996 ) . We investigated the effect of alter- 
ing the slope and found that a change in slope of ~ 0.2 does 
not significantly alter our conclusions. 

The conversion from spatial to angular cross-correlation 
function is obviously strongly dependent on the radial selec- 
tion functions, 4>(x), of the two populations. For the QSOs 
we simply use the redshift distribution of the QSOs used in 
our analysis above (Section ^|), that is, we integrate Lim- 
ber's equation over all the QSOs in the sample. To model 
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Figure 6. Results for the QSO cross-correlation with colour selected galaxies. The angular cross-correlation function for a) red and b) 
blue galaxies. The clustering amplitude within 120" for individual QSOs with c) red and d) blue galaxies as a function of redshift, and 
the redshift binned clustering amplitude for e) red and f) blue galaxies. 
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Figure 7. The model redshift distribution from Metcalfe et al. 
(1996) for a limiting magnitude of 6j = 23. 



Table 3. Predicted clustering amplitudes, u>(9 < 120") and A qg 
(assuming the given value of 7), for 2 < 1.5 QSOs and foj < 23 
galaxies in a number of different models, a gives the significance 
of the difference between the models and the observed clustering 
for all QSOs with 2 < 1.5. No correction is made for the integral 
constraint. 



r 7 e w(0 < 120") A qg a 
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1.8 


-1.2 


0.0430 
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the galaxy selection function, we use the N(z) derived from 
the galaxy evolution model of Metcalfe et al. (1996) which 
is consistent with the redshift distribution of galaxies in 
t he deep Keck Tele scope galaxy redshift survey to B — 24 
(Cowie et al. 1995). These models contain an exponential 



rise in star formation with look-back time and use the evo- 
lutionary tracks of Bruzual & Chariot (1993). The N(z) dis- 
tribution from this model is shown in Fig. It should be 
noted that the luminosity evolution produces a considerable 
high redshift tail to the distribution, as found in the sample 
of Cowie et al., where significant numbers of galaxies were 
found above 2 = 1. 

Therefore, the only free parameters we have are the 
clustering correlation length ro and the evolutionary param- 
eter e, which we attempt to constrain from the QSO-galaxy 
cross-correlation results. 



4.2 A comparison of models and data 

We calculate the expected clustering amplitude, u(6 < 
120"), for a number of different cases which are shown in 
Table [| We first use the known present day galaxy-galaxy 
auto-correlation function with ro = 6 ft" 1 Mpc and 7 = 1.8. 
With clustering evolution which is stable in comoving co- 
ordinates this model is inconsistent with the data for all 
2 < 1.5 QSOs at the 3.5a level. Clustering which is stable 
in proper coordinates is rejected at the 2.6a level and clus- 
tering evolution which is consistent with linear theory is re- 
jected at the 2.2a level. If we allow for an integral constraint 
of ~ 0.006, which of course assumes a clustering amplitude 
similar to the auto-correlation amplitude of 6.7 < 23 galaxies, 
the significance of the results are only changed by ~ 0.3<r. 
This correction would be considerably smaller if we used 
the measured cross-correlation amplitude to estimate an in- 
tegral constraint. In fact the integral constraint would be 
negative, increasing the significance of the above rejections. 
Results from three other sets of parameters are also shown in 
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Figure 8. u)(9) calculated from the models discussed in the 
text compared to the results from all QSOs with 2 < 1.5 cor- 
related with fej < 23 galaxies. The solid curves show models 
with comoving clustering evolution (e = —1.2) and amplitudes 
of ro = 14 h^ 1 Mpc (upper) and 6 h~ x Mpc (lower). The dot- 
ted curves show models with linear theory clustering evolution 
(e = 0.8) and amplitudes of ro = 14 h~ 1 Mpc (upper) and 
6 h^ 1 Mpc (lower). 
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Figure 9. a) w(0 < 120") and b) w(0 < 240") calculated from 
the models discussed in the text compared to the results from 
all QSOs as a function of redshift. The solid curves show models 
with comoving clustering evolution (e = —1.2) and amplitudes of 
ro = 14 h~ 1 Mpc (upper) and 6 h _1 Mpc (lower). The dotted 
curves show models with linear theory clustering evolution (e = 
0.8) and amplitudes of ro = 14 h _1 Mpc (upper) and 6 h^ 1 Mpc 
(lower) . 



Table ^. With a clustering amplitude similar to that of the 
cluster-cluster auto-corrclation fun ction ro = 14 ft -1 Mpc 
and 7 = 1.8 (Dalton et al. 1994), the lowest amplitude 
found (assuming linear theory growth) is 5a too high com- 
pared with the data. We then try a model which fits the 
galaxy-cluster cross-correlation functio n, ro = 8.8 h~ x Mpc 
and 7 = 2.2, flLilje fc Efstathiou 198S| ). This model, which 
might be a better estimate of our expected signal if QSOs 
were found only in rich clusters, is comfortably rejected for 
any reasonable rate of evolution. Finally we test a model 
with very low clustering amplitude, ro = 2 h^ 1 Mpc, simi- 
lar to that found in the Canada- France Redshift Survey ( [Ee] 
Fevre et al. 1996|) for 7ab < 22.5 galaxies with a median 



redshift of 0.56. This is clearly a better fit to our results, 
being discrepant by only ~ 1.5cr. Ideally we would attempt 
to find the best fit values for ro and e using maximum likeli- 
hood techniques. However, as the measured cross-clustering 
amplitude is marginally negative, these best fit values would 
have no physical meaning. Fig. |^ demonstrates that on scales 
> 120" the data are also inconsistent with the models, out 
to ~ 8'. The cross-correlation is, on average, negative up to 
this scale and inclusion of an integral constraint of the mag- 
nitude discussed above (Section Eh has no significant effect 
on this statement. 

The models are also plotted as a function of redshift 
in Fig. § Fig. |a is for u(0 < 120") while Fig. |b is for 
u>(9 < 240"). These first demonstrate the relative stability 
of using a 120" radius aperture to calculate the clustering 
amplitude. Secondly, they show the redshift at which the 
models can rejected. In Fig. ^b, all the models plotted are 
inconsistent with the data for redshifts less than z ~ 1.5. 



4.3 The effect of gravitational lensing 

Our results suggest that there could be a lensing compo- 
nent to the cross-correlation of faint, 6j < 23, galaxies and 
optically selected QSOs, as a small anti-correlation was mea- 
sured between the two populations. Similar anti-correlations 
between brighter galaxies and UVX QSO candidates have 
been interpreted in terms of gravitational lensing by Croom 
(1997). This anti-correlation is due to the two competing 
effects of magnification and area distortion. If the number- 
magnitude relation of the lensed objects is flat, there are 
only a small fraction of objects that are added to the sam- 
ple through amplification, so that the area distortion wins, 
reducing the perceived number density. 

Here we attempt to model the effect of the lensing con- 
tribution on this result. Using the redshift distributions of 
the galaxy and QSO populations we construct Monte-Carlo 
simulations to calculate the expected correlation signal as 
a function of redshift. These simulations model the individ- 
ual galaxy potentials as isothermal spheres, and take into 
account the clustering of galaxies by also adding a constant 
density plane. Wu et al. (1996) calculate that the maximum 
matter surface density contributed by large-scale structure 
is ~ 0.01 — 0.02/i g cm -2 assuming that the amount of mat- 
ter in galaxies is fi g ~ 1. Typical velocity dispersions for 
galaxies are ~ 200kms~ 1 , so we produce Monte-Carlo sim- 
ulations with velocity dispersions of 200 and 400kms _1 and 
constant density planes with E = 0.0 and 0.02/i g cm -2 . We 
use a value of 0.28 for the slope of the QSO integral number- 
count slope. Fig. [h] shows how the clustering amplitude in 
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these simulations changes with redshift. The extreme case 
of E = 0.02/i g cm" 2 and a = 400kms~ 1 , which can con- 
tribute ~ —0.02 to the correlation amplitude, can account 
for most of the negative signal found for QSOs at z < 1.5, 
while it cannot account for the negative correlation found 
at higher z. Fig. [l0^ 3 shows the u){6) predicted by the lens- 
ing models, for a combination of all the QSOs at z < 1.5, 
compared to the data. Correcting the data to allow for this 
lensing effect would not give the correlation function a posi- 
tive amplitude. Lensing has no effect on the rejection of high 
clustering amplitude models, but does make the models with 
low clustering amplitude more consistent with the data. 



5 DISCUSSION 

It is clear that there is strong disagreement between the clus- 
tering results presented here and a range of models. This is 
demonstrated by Fig. |^. At z < 1.5 the correlation within 
120" is u(9 < 120") = -0.027 ± 0.020, and at z > 1.5, the 
combined X-ray and optically selected sample has a clus- 
tering amplitude of u(6 < 120") = -0.028 ± 0.021, almost 
identical to that found at lower redshift. If we assume that 
a non-zero cross-correlation is only produced by physically 
associated QSOs and galaxies then the above results give us 
strong constraints on the environments of radio-quiet QSOs 
up to z ~ 1 — 1.5. Out to this redshift we still expect to see 
significant numbers of galaxies with a limiting magnitude 
of bj = 23. The local luminosity function has Mj ~ —19.7 
(e.g. Ratcliffe 1996) which is equivalent to 6j = 23 and 24 for 
z — 1.0 and z — 1.5 without K-corrections. A combination of 



K-correctio n and reasonable luminosity evolution (Pozzetti 
et al. 1996 ) leaves these values largely unchanged (partic- 
ularly for late type galaxies) and thus at z — 1 we reach 
~ M* and at z — 1.5 about 1 mag brighter than M* . We 
can firmly reject the hypothesis that the QSOs in our data 
exist in rich environments similar to galaxy clusters. This is 
true for a wide range of evolutionary rates (parameterized 
by e). A scenario in which QSOs have poorer environments, 
similar to 'normal' galaxies, is a more acceptable alterna- 
tive, although this still gives a clustering amplitude which 
is inconsistent with the data at the 2 — 3<r level. 

We see a significant anti-correlation in our data, which 
cannot be explained by the conventional clustering models 
used above, but the data give us various clues as to the 
source of this signal. First, we find that the optically se- 
lected QSOs show more anti-correlation than the X-ray se- 
lected samples, particularly at high redshift. Secondly, when 
we sub-divide the galaxy samples on the basis of colour 
(for only two fields), we find that the red galaxies show a 
stronger anti-correlation than the blue galaxies. It has been 
established that the auto-correlation amplitude for faint red 
galaxies is significa ntly higher than for faint blue galaxies 
(Roche et al. 1996). This is easily understood if the redder 



galaxies lie predominantly at low redshift while the high red- 
shift tail of the galaxy distribution is composed of the blue 
population. Of course, early type galaxies are also more nu- 
merous in rich structures such as galaxy clusters. It therefore 
appears that QSOs tend to show an anti-correlation with 
galaxies in rich environments at low redshift. 

The possible causes of this signal are either obscura- 
tion by high densities of dust in rich structures or gravi- 
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Figure 10. Estimated angular cross-correlation function of QSOs 
with faint, bj < 23, galaxies due to lensing. a) lu(8 < 120") as 
a function of QSO redshift. The correlation amplitude is mea- 
sured within 120" from Monte-Carlo simulations, with velocity 
dispersions 200 and 400kms~ 1 and planes with S = 0.0 and 
0.02/i g cm _2 mass surface densities, b) w(9) for the four differ- 
ent lensing models for all QSOs with z < 1.5 compared to the 
measured u(9). 
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tational lensing which can cause an anti-correlation if the 
QSO number-counts slope is flat (e.g. Narayan 1989). The 
optical QSOs in the sample are in the magnitude range 
where the number counts are flat, ~ 0.3, compared to the 
critical slope of 0.4 which is the slope required for zero- 
correlation. The X-ray selected QSOs, in contrast, have a 
logiV — logS 1 slope of 2.52 zfc 0.3 at bright magni tudes turning 
over to 1.1 ± 0.9 ( Gcorgantopoulos et al. 1996] ). The critical 
logiV — logS slope is 1.0 so we expect no anti-correlation 
(although it should be noted that the errors in the slope at 
faint fluxes are large). Above we show that even extreme esti- 
mates of galaxy masses could only produce anti-correlations 
of u>(8) ~ —0.02 to —0.03. This extreme model can just 
account for the anti-correlation measured. Lensing could re- 
duce the inconsistency between our results and the lowest 
amplitude clustering models to only ~ la, while the in- 
consistency with the strongly clustered models is still very 
significant. Lensing does not appear to explain the larger 
anti-correlation found at z > 1.8 between the bj < 23 galax- 
ies and optical QSOs, although it should be noted that this 
signal is due to a relatively small number of QSOs. A final 
possibility is that QSOs at low redshift might avoid rich en- 
vironments. However, we have only a small number of QSOs 
at low redshift which are in fields that contain colour infor- 
mation (7 QSOs a t z < 1), and also, other observations 
(Smith et al. 1995) find a low but positive correlation be- 
tween QSOs and galaxies at low redshift. 

Finally we note that deep wide field CCD imaging will 
allow this study to be developed further. Deep, B < 26, 
imaging in an area of ~ 2 sq. deg. will be presented in the 
second paper in this series (Croom & Shanks 1998 in prepa- 
ration). This data will allow us to probe QSO environments 
at higher redshift z ~ 2, and determine if the present study 
may have missed small, faint galaxies associated with lower 
redshift QSOs. 



6 CONCLUSIONS 

We have carried out an investigation of the environments 
of a large sample of QSOs (~ 150), using deep AAT pho- 
tographic plates in 5 independent areas. The QSOs were 
selected by both optical colour techniques and deep X-ray 
exposures from ROSAT. We draw the following conclusions: 

1. Most QSOs show zero or negative correlations. For a 
combination of all QSOs at z < 1.5 we find u(6 < 120") = 
-0.027 ± 0.020. Only a small number of QSOs (~ 4 - 5) 
show significant positive correlation with 6j < 23 galaxies; 

2. Optically and X-ray selected QSOs show marginally 
different properties, with the optical sample being more anti- 
correlated than the X-ray sample; for z < 1.5 the X-ray 
selected sample gives u(6 < 120") = -0.018 ± 0.022, while 
the optically selected sample gives u)(9 < 120") = —0.051 ± 
0.027. At high redshift, z > 2, the optical QSOs become 
more anti-correlated with galaxies. 

3. Analysis of the two fields with colour information sug- 
gests that the anti-correlation is stronger in the red galaxy 
population than in the blue population. 

4. The cross-correlation results are compared with mod- 
els which include clustering evolution. Models which have 
clustering amplitudes similar to that of galaxy clusters (ro ~ 
14 h' 1 Mpc) with a wide range of evolutionary rates are 



ruled out. The data are more consistent with a low cluster- 
ing amplitude, similar to that of the galaxy auto-correlation 
function, although ro = 6 hT 1 Mpc combined with linear 
theory evolution still gives too high an amplitude at the 
2a level. Models with ro = 2 h^ 1 Mpc are only ~ 1.5a 
above the observed cross-correlation. If gravitational lens- 
ing causes some of the measured anti-correlation, the dis- 
agreement with the low amplitude clustering model could 
be reduced to ~ la. 

5. We suggest that the anti-correlation found could be 
due to gravitational lensing, particularly as the optical QSOs 
and red galaxies show the strongest anti-correlation. An ex- 
treme lensing model could account for most of the measured 
anti-correlation at z < 1.5. However, even this extreme lens- 
ing model does not significantly affect our conclusions when 
comparing the data to clustering models. 

The results presented here clearly reject the notion that 
radio-quiet QSOs exist in rich cluster-like environments up 
to z ~ 1 — 1.5. Beyond this redshift we require deeper images 
to determine the environments of QSOs, these observations 
will be important in the interpretation of new QSO large- 
scale structure measurements from the next generation of 
QSO surveys (e.g. the 2dF QSO Redshift Survey, Smith et 
al. 1997). 

Our results are consistent with those of Boyle & Couch 
(1993) who find no correlation between R ~ 23 galaxies and 
QSOs in the redshift range 0.9 < z < 1.5. Ellingson et al., 
(1991) find a marginal positive signal at z < 0.6, which is 
consistent with our data (although we have a small number 
of QSOs at these low redshifts). Our data appear to be in- 
consistent with those of Hutchings et al., (1995) who find 
that a sample of 9 out of 14 QSOs are associated with com- 
pact groups of star forming galaxies, although these objects 
lie in a QSO super-cluster and so may be a unrepresenta- 
tive sample existing in a particularly rich environment. We 
conclude that our results continue to uphold the hypothesis 
that QSOs are not strongly biased with respect to the galaxy 
population, but are more likely to trace the distribution of 
normal galaxies. 
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